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We consider braneworld scenarios including the leading correction to the Einstein-Hilbert action 
suggested by superstring theory, the Gauss-Bonnet term. We obtain and study the complete set 
of equations governing the cosmological dynamics. We find they have the same form as those in 
Randall-Sundrum scenarios but with time-varying four-dimensional gravitational and cosmological 
constants. Studying the bulk geometry we show that this variation is produced by bulk curvature 
terms parameterized by the mass of a black hole. Finally, we show there is a coupling between these 
curvature terms and matter that can be relevant for early universe cosmology. 



In recent decades developments in cosmology have 
been strongly influenced by high-energy physics. A re- 
markable example of this is the inflationary scenario and 
all its variants. This influence has been growing and 
becoming more and more important. Today it comes 
from developments in string and M theories and the new 
scenarios they are providing for cosmology. In particu- 
lar, the study of scenarios where the spacetime has non- 
compact extra-dimensions has produced intense activity. 

The most popular model of that class is the one pro- 
posed by Randall and Sundrum §,| (RS), motivated 
by orbifold compactification of higher dimensional string 
theories, in particular, by the dimensional reduction of 
eleven-dimensional supergravity in R 10 x S 1 /I2 intro- 
duced by Hofava and Witten Q. The picture coming 
from the RS model is one with all matter and gauge 
fields, except gravity, confined in a three-brane embedded 
in a five-dimensional (5D) spacetime with Z2 symmetry, 
and where the zero-mode of the Kaluza-Klein (KK) di- 
mensional reduction is localized, reproducing Newtonian 
gravity in the weak field approximation. 

On the other hand, it is a general belief that Einstein 
gravity is a low-energy limit of a quantum theory of grav- 
ity which is still unknown. Among promising candidates 
we have string theory which suggests, that in order to 
have a ghost-free action, quadratic curvature corrections 
to the Einstein-Hilbert action must be proportional to 
the Gauss-Bonnet term This term also plays a fun- 
damental role in Chern-Simons gravitational theories [^). 
However, although being a string-motivated scenario, the 
RS model and its generalizations |6) do not include these 
terms. From a geometric point of view, the combination 
of the Einstein-Hilbert and Gauss-Bonnet term consti- 
tutes, for 5D spacetimes, the most general Lagrangian 
producing second-order field equations |?J (see also Q). 

These facts provide a strong motivation for the study 
of braneworld theories including a Gauss-Bonnet term. 
Recent investigations on this issue have shown Q that 
the metric for a vacuum 3-brane (domain wall) is, up 



to a redefinition of constants, the warp-factor metric of 
the RS scenarios. The existence of a KK zero-mode lo- 
calized on the 3-brane has also been demonstrated (h]] 
(see also [|ll]jlj]). Properties of black hole solutions in 
AdS spacetimes have been studied in The cos- 

mological consequences of these scenarios are less well 
understood. This issue has been studied in [l5Ul^| in the 
case of a 2-brane model Q and in |[7| for a single-brane 
model However, in both cases only simple ansatze 
for the 5D metric (written in Gaussian coordinates as 
in |l8|]) were considered, e.g. the separability of the met- 
ric components in the time and extra-dimension coordi- 
nates. One can see that this assumption is too strong 
even in RS cosmological scenarios p8| , where they lead 
to a very restrictive class of cosmological models, not 
representative of the true dynamics. Other works with 
higher-curvature terms in braneworld scenarios are con- 
sidered in pS| . 

In this paper we obtain the equations governing the 
dynamics of Friedmann-Robertson- Walker (FRW) cos- 
mological models in braneworld theories with a Gauss- 
Bonnet term. This includes the derivation of the appro- 
priate junction conditions for this kind of theory. Then, 
we investigate the new dynamical cosmological behaviour 
and how it is related to the geometry of the bulk. 

Our starting point is the following action 

+ J d 4 x^j (C m - 2A) , (1) 

where Ceh = R — 2A is the Einstein-Hilbert Lagrangian 
with a negative cosmological constant, A < 0, Cgb is the 
Gauss-Bonnet correction 

C GB =R 2 ~ AR ab R ab + R abcd R A bcd , (2) 

and matter fields (C m ) are confined to a 3-brane with 
Z2 symmetry. Moreover, A is a constant that coincides 
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with the brane tension in the limit C m = = a. Ob- 
jects corresponding to the brane are written with a tilde 
to be distinguished from 5D objects; g is the metric de- 
terminant and R, Rab, and Rabcd are the scalar cur- 
vature and Ricci and Riemann tensors respectively; k is 
the 5D gravitational constant and physical units in which 
c = 1 are assumed. The sign of the fundamental constant 
a must be positive according to the expansions carried 
out in string theory . If we write the unit normal to 
the 3-brane as n = dw (Gaussian coordinates), the field 
equations are 

Gab + kg A B + aH A B = n 2 h^fJAB + T A b] S(w) , (3) 

where Uab is the second-order Lovelock tensor 

Hab = RRab — 2Ra c Rbc — 2.R D Racbd 

+Ra CDE Rbcde — j9ab£gb , 

Tab is the energy-momentum tensor describing the mat- 
ter confined on the 3-brane (TabtiP = = gABn B ). 
Since we are interested in the cosmological consequences 
of the theory we take Tab to be of the perfect-fluid type 



Tab = (p + p)uau b + P9ab , 
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with u A , p, and p, being the fluid velocity (u a ua 
energy density and pressure respectively. 

Homogeneous and isotropic cosmological models can 
be described by the 5D line element |TJ| (i,j,... = l, 2, 3) 

ds 2 = -n 2 (t, y)dt 2 + a 2 (t, y)h i3 dx i dx j + b 2 (t, y)dy 2 , (4) 

where y is the fifth dimension coordinate (the brane 
is located at y = 0) and hij is a 3-dimensional maxi- 
mally symmetric metric for the surfaces {t = const., y = 
const.}, whose spatial curvature is parameterized by 
k = — 1, 0, 1 . Every hypersurface y — const, has the met- 
ric of a FRW cosmological model. 

The first step to solve the field equations (||) is to study 
them in the bulk (y ^ 0). Following [Q we have found 
that a set of functions {a(t,y),b(t,y),n(t,y)} constitute 
a solution of the field equations in the bulk provided the 
following equations are satisfied f2G] 
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where C is an integration constant and 
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*(*»!/) = -a — - S— + — • 
n z a z b z a z a z 

The simplicity of equation (j^ is remarkable and so is the 
fact that the Einstcin-Hilbert and Gauss-Bonnet terms 



depend on the same argument, namely <!> . It would be 
interesting to know whether or not this happens also for 
higher-order terms in D > 5 spacetimes, i.e., whether the 
dependence of the nth-order Lovelock term is of the form 

To find the equations for the 3-brane we need to study 
the junction conditions of our theory, which will be dif- 
ferent from those in Einstein gravity |2lJ . We know that 
the 5D metric must be continuous across the brane and 
that there will be jumps in its normal derivatives due 
to the energy-momentum distribution on the brane [see 
Eq. (H)]. The Z2 symmetry makes the metric invariant 
under the transformation y — > —y and hence any metric 
component can be written as A(t,y) = A(t, \y\). There- 
fore, the normal derivatives are given by 

A'(t, y) = d ]y] A(t, M)sign(y) = A x (t, |y|)sign(y) , 

A"(t,y) = d lyl A 1 (t,\y\)+2A 1 (t,\y\)S(y). 

The first derivative has a finite discontinuity across 
the brane but its square is continuous by virtue of 
the Z2 symmetry. The second derivative has a non- 
distributional part, the first term, and a distributional 
part whose coefficient gives the jump across the brane, 
namely [A'] = A'(t,0 + ) - A'(t,Q-) = 2A 1 {t,Q) . The 
value of this jump can be obtained by integrating the 
tt-component of the field equations (||) across the brane, 
that is (5(w) = b- 1 8(y)) 



lim 



dy (G tt + A 5tt + aH tt ) = n 2 (p + , 



where from now on the subscript "o" will denote the value 
of the corresponding quantity on the brane. Then, the 
equation for the jump 5i(t, 0) is 
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This is a cubic equation for the discontinuity a\ which, 
for sufficiently small a, has only one real solution, the 
other two being complex. Therefore, if we require our 
cosmological equations to have the right a — > limit we 
are left with only one solution. Introducing it into equa- 
tion (Q) will lead to the Friedmann equation for the brane. 
Equation (^|) is a quadratic equation in the square of the 
Hubble function H = a /(n a ), but only one solution 
has the correct limit a — > 0. Remarkably, the Friedmann 
equation we get from that solution can be written in the 
same form as in the RS cosmological scenario (compare 
with the generalized Friedmann equations in |2^j23|) 
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where the 4D gravitational coupling and cosmological 
constants, k* and A*, are now time dependent. The 
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only difference with the corresponding RS equation is 
that here we do not have explicitly the dark radiation 
term ||| proportional to a~ 4 . It is actually included, in 
a non-linear way, in A* [see Eq. (|[) below]. The time- 
dependent fundamental "constants" k 2 and A* change 
in time as functions only of the scale factor a . Their 
explicit form is 
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(8) 



where we have introduced cr, the function containing 
the dependence on a a . Moreover, k 2 is the 4D gravi- 
tational coupling constant, the one that appears in the 
computation of Newton's law pT| . For small a we re- 
cover the Friedmann equation in RS braneworlds (see, 
e.g., J2^j23|). Actually, at zero order in a, A* = A/jg + 
3Ca~ + 0(a) , recovering the dark radiation term Ca~ 
in the Friedmann equation of RS brane- worlds. 

The dependence of k 2 and A* on a has been plotted 
in Fig. [l]. There are two possible behaviours according to 
whether A* has a minimum [Fig. 0(a)] or not [Fig. 0(b)] . 
In both cases A* tends to infinity for small a a and to a 
constant 



Aoo = 



2,i 



(1 



1 



aA) 



when a a — > oo. Then, A* tends to zero for the following 
critical value of A 



A, = 
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a k 2 



(1-Vl + foA). 
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corresponds to a 



Minkowskian brane jl],|| (see J9|). From Eq. (|) we must 
have 1 + 4aer > 0, hence from the behaviour for big and 
small a we deduce that a|A| < 3/2 and aC > 0. On the 
other hand, the behaviour of k 2 : when C ^ 0, is the same 
independently of the value of the parameters; it tends to 
zero for small a and to k 2 [see Eq. (^)] when a a — > oo . 

As in standard cosmology, the Friedmann equation (0) 
together with the energy-momentum tensor conservation 
equations [a consequence of the divergence-free character 
of the left-hand side of (0)], 

p=-3(p + p)H, 

and a barotropic equation of state p = p(p), describe 
completely the cosmological dynamics on the brane. 

At this point we do not know much about the ge- 
ometry of the bulk. For perfect-fluid FRW cosmolog- 
ical models in RS braneworlds, the bulk turns out to 
be the Schwarzschild-AdS 5D spacetime (see p4] , p5| ). 
Since in our case we have different gravitational equa- 
tions [Eqs. ([| ||)] we do not expect the same bulk. The 
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FIG. 1. Variation of the effective 4D gravitational and cos- 
mological constants, k 2 t and A*, with respect to the scale fac- 
tor a [Eq. (@)]. (a) A* has a minimum, (b) A* is a monotonic 
function of a a . 



knowledge of the bulk geometry is important to under- 
stand the physical meaning of the integration constant 
C appearing in the Friedmann equation (f^,^|). Assum- 
ing that the fifth dimension is static |lS| ], 6 = 0, we have 
found that there is a coordinate change, analogous to 
the one used in RS scenarios (26|, that brings the line 
element (0) to the following static form 



ds 2 = -f(R)dT 2 



dR 2 

f(R) 



+ R 2 (d X 2 + Z 2 k (x)dn 2 ), (9) 



and dfl?. 



where E_i = sinhx , So = x > an d Si = sinx 
is the unit 2-sphere metric. In RS scenarios, Kraus |24|] 
and Ida p5[ , found the Schwarzschild-AdSs bulk. When 
a Gauss-Bonnet is present, the solution of the field equa- 
tions (j|) for the line-element (||) is given by (see [Q for 
the k = case and [(Til for any k) 



f(R) 




(10) 



where 14 is the volume of the 3-surface {T — const., R = 
1}. In this picture, the 3-brane is a hypersurface that 
expands or contracts according to a law R = R(T) , where 
R coincides with the scale factor, R(T) — a(t,y), and t 
is a proper time (n = 1). More important, the constant 
C is related to the black hole mass by the relation 



C = 



k M 
~3Vk 



Therefore, the string theory prediction that a must be 
positive implies, using aC > 0, that the mass of the 
black hole must be positive. Conversely, if we require the 
avoidance of a naked singularity, M > H], then the 
constant a must be positive. 

As one would expect, if the mass of the black hole tends 
to zero the bulk Weyl tensor vanishes. This means that 
the effective time- variation of k 2 is due to the coupling of 
bulk Weyl curvature terms with matter. From a geomet- 
rical point of view, it is a consequence of the structure 
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of the curvature quadratic terms in the Gauss-Bonnet 
Lagrangian. 

Let us now study what is the behaviour of FRW mod- 
els in this theory. We start with the most simple case, 
when there is no black hole present in the bulk, which 
means M = 0. In this case the dynamics is completely 
equivalent to that of RS braneworlds, widely discussed in 
the literature (see [^3],^7| for an exhaustive study), where 
fi and Aqo are the 4D gravitational and cosmological con- 
stants respectively. Moreover, one can check that in this 
case the conditions for inflation are the same as those 
described in p8) . 

To study the case M ^ we assume a linear equation 
of state, p = (7 — l)p, which implies p oc a~ 37 . For late 
times, a a 3> 1, and as in RS braneworlds, the behaviour 
is the same as in standard cosmology, a ~ t 2 ^ 3l \ but 
with the modified 4D constants k and A^. This is ex- 
pected since it corresponds to the low-energy regime. At 
high energies, or early times (a Q <C 1), things are very 
different. In the case 7 > 1 we find that H 2 as a^ -67 , 



hence, a a 

n 



£1/(37-2) , The generic behaviour for 7 < 1 
is H 2 as a~ 2 , therefore a a as t. Then, when 7 < 1 the 
dynamics is independent of the particular equation of 
state. This results show that the dynamical behaviour 
in our theory is completely different from the general 
relativistic case (a as i 2 /( 37 )), except for the radiation 
case (7 = 4/3), when they remarkably coincide. It 
is also different from the behavior in RS braneworlds 
(o„ as tV(37)). 

In conclusion, we have studied the cosmological dy- 
namics in a theory that generalizes the RS scenario [jl]]^] 
by taking into account the Gauss-Bonnet higher-order 
curvature term [|]J^]. Using this fact, we have shown 
that the cosmological equations can be seen as those of 
RS scenarios but with time-dependent 4D gravitational 
and cosmological constants. Studying the 5D geometry 
of our model we have found that the time variation of the 
constants is parametrized only by the mass of a black hole 
in the bulk. In the case of the gravitational constant, the 
time dependence is introduced through the coupling be- 
tween bulk curvature terms and matter. Finally, we have 
shown how the higher-order curvature terms present in 
our theory are dominant at high energies and change the 
cosmological dynamics at early times. Hence, they can 
provide alternative cosmological scenarios for the study 
of unsolved cosmological problems. 
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